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We propose to use the complex-range Gaussian basis functions, {r'e"'-^^*'^-"-'^'^'""-' Yim(r); 
r„ in a geometric progression}, in the calculation of three-body resonances with the complex- 
scaling method (CSM) in which use is often made of the real-range Gaussian basis functions, 

^ , {r' e-('-/'-") F,„.(f)}, that are suitable for describing the short-distance structure and the 

. asymptotic decaying behavior of few-body systems. The former basis set is more powerful 

than the latter when describing the resonant and nonresonant continuum states with highly 
oscillating amplitude at large scaling angles 6. We applied the new basis functions to the CSM 
QQ , calculation of the 3a resonances with J = 0^, 2^ and 4^ in ^^C. The eigenvalue distribution 

of the complex scaled Hamiltonian becomes more precise and the maximum scaling angle 
becomes drastically larger (6'max = 16° — >■ 36°) than those given by the use of the real-range 
Gaussians. Owing to these advantages, we were able to confirm the prediction by Kurokawa 
and Kato [Phys. Rev. C71, 021301 (2005)] on the appearance of the new broad Oj state; 
we show it as an explicit resonance pole isolated from the 3q continuum. 

^ . §1. Introduction 

C . 

The complex scaling method (CSMjIl'''21'''3''31''Eli is a very powerful tool to in- 
vestigate many-body resonances that recently attract much interest in experimental 
and theoretical developments in unstable nuclear physics where bound, resonant and 
l/^ ' continuum states are energetically degenerate and must be treated simultaneously. 

CnI ; Using the CSM, one can directly obtain the energy Ey. and the decay width F of 

a many-body resonance by solving the eigenvalue problem for the complex scaled 
' Schrodinger equation with a scaling angle 0, [H{0) — E{9)\F{6) = 0, where the 

. boundary condition of the many-body outgoing wave is automatically satisfied for 

the resonance, giving E = E^ — iF/2, that is in principle independent of 6. A 
merit of the CSM is that the resonant eigenstates are described by square-integrable 
wave function. The CSM and its application to the nuclear physics problems are 
extensively reviewed in Ref. ^ and references therein. 

In the CSM, there is a limitation of the scaling angle 9 due to the analyticity of 
the Hamiltonian. Furthermore, in practical calculation, one often meets a difficulty 
in solving resonant states with large decay width since the complex scaled Hamil- 
tonian is diagonalized with a limited number of basis functions. A set of real-range 
Gaussians, {r' e"'-''/'""-' lim(r); ^'n in a geometric progression}, are often 
employed as the basis functions. 

Thus, one of the purposes of the present work is to propose the use of the 
complex-range Gaussian basis functions, {r' e"*^^^*'^)*^''/''") Yim(r)}^^ in the CSM 
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calculation of three-body resonances so as to overcome the above difficulty for the 
large 6. Owing to the oscillating component, the space of the new function set is 
much larger and the linear dependence between basis members is much smaller than 
those in the basis set of real-range Gaussians. This improves the quality of the CSM 
calculation significantly and increases the possible scaling angle drastically. 

One of the most intensively studied nuclei using the three-body CSM is ^^C 
nucleus as the 3a-cluster system. The CSM has especially been useful to investigate 
the 3a resonance structure in its excited states since the method satisfies the cor- 
rect boundary condition for the three particles decaying into three-body continuum 
states. Such CSM studies of are reviewed in Refs. [TT]) and \T2^ and references 
therein. Among the studies, Kurokawa and KatcP^'l^ succeeded in thoroughly cal- 
culating the energies and decay widths of the 3a resonant states in ^^C with the 
total angular momentum J = to 5. 

Table I. New results of the recent calculatiorP^''^ and experimenlP^ on the 0;^, and 2^ 
states in ^^C. The excitation energies (Ex) and decay widths (F) are given in MeV. See also 
the footnote below. 



12fj 
J" 




EXF 


131 


Ex 


r 


Ex 


r 




8.95 


1.48 


9.04 ± 0.09 


1.45 ±0.18 


ot 


11.87 


1.1 


10.56 ±0.06 


1.42 ±0.08 


2t 


9.57 


1.1 


9.84 ± 0.06 


1.01 ±0.15 



In Table I, it is interesting to see the recently accomplished reasonable agree- 
ment between the results by the calculationP^'l^ and the observatioEpJ on the Og", 

and 22 states. Especially, Kurokawa-Kato's prediction of the new Oj^ state hav- 
ing a large width is of importance to understand the new experimental results in 
Ref. [EDO Their CSM calculation was performed at the scaling angle of 16° that 
is the largest angle available in the calculation using the real-range Gaussian basis. 
Since the angle is not enough for separating the low-lying broad Oj^ resonance from 
the continuum eigenvalues, they made an extrapolation by applying the method of 
analytic continuation of the coupling constant (ACCC^(^''^'I^ combined with the 
CSM (ACCC-I-CSM) in order to derive the complex energy of the Oj^ resonance. 
But, AraP3 reported that, the 0^ resonance obtained in Ref. Ilip is missing from his 
calculation based on the microscopic i?-matrix method for the ®Be(0+, 2"'", 4+) + a 
two-body scattering problem. However, in the calculation of Ref. [T7|) . the ^Be was 
described with the bound-state approximation employing only four different tem- 
pered Gaussian functions for the a-a relative motion. 

*-* As for the O""" data in Fig. 8(a) in Ref. I13[) . we employ the interpretation by the authors 
that there are two 0"'' resonance peaks as summerized in Table I and that the resonances may 
correspond respectively to the 03" and Oj states described in Ref. Ilip . The authors showed another 
interpretation to regard the peaks as a single peak at Ex = 9.93 ± 0.03 MeV with a width of 
2.71 ±0.08 MeV. 
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Thus, the second purpose of the present paper is to apply the CSM with the 
complex-range Gaussian basis functions to the 3a resonances in ^^C and examine 
the results in Refs. Ill|) and IT^ . Since the function space of the present basis set 
is very large, we have the following advantages: i) distribution of the eigenvalues of 
the complex scaled Hamiltonian becomes much more precise than those obtained in 
the literature, and ii) the scaling angle is drastically increased from 9 = 16° up to 
36° that is enough large to separate explicitly the Oj^ resonance pole from the 3a 
continuum eigenvalues. 

The present paper is organized as follows: in §2, we introduce the complex-range 
Gaussian basis functions and incorporate them into the framework of CSM. In §3 
we apply the method to the 3a resonances in ^^C and compare the result with that 
obtained in Refs. TP) and 1 1 2,p . Summary is given in §4. 

§2. Method 

2.1. Three-body complex scaling method 

In many cases of the CSM studies, the three-body wave function is expanded in 
terms of the real-range Gaussian basis set with ranges in a geometric progression. In 
this work, we propose to use the complex-range Gaussians in the three-body CSM 
calculations. 

We explain it, as an example, taking the case of ^^C(= a-|-a-|-a) on the basis of 
the orthogonality condition model (OCM)^^' for the 3a system. The extension from 
the real-range Gaussian to the complex-range ones in other three-body systems is 
straightforward. We take all the three sets of Jacobi coordinates (Fig. 1), ri = X2— X3 
and Ri = xi — ^(x2 + X3) and cyclically for (r2,R2) and (r3, R3), Xj being the 
position vector of ith particle. 




i=1 i=2 i=3 

Fig. 1. Three sets of the Jacobi coordinates for three a particles. 

The Hamiltonian is written as 

3 3 

H = Y,ti-TG+Yl ^c,a{n) + F3«(ri, r2, r3) + Fpauli- (2.1) 

1=1 i=l 

The operators and Tq stand for the kinetic energies of ith a particle and the 
center of mass motion, respectively. Vaa is the a-a potential and l^auii is a pseudo 



potential representing the Pauli principle between a clusters. The 3a potential Y^a 
is introduced if necessary. These potentials is explained in §3.1. 
In the CSM, the radial coordinates are transformed by 

n^r,e'\ R,^Rie'^. (2.2) 

The transformed Hamiltonian is denoted by H{9). We solve the equation 

[H{e) - E{e)]^{e) = o (2.3) 

by expanding •P'(6') in terms of the totally symmetric L^-integrable three-body basis 
functions {^j, 7 = 1, • • • , 7max} : 

O'max 

m = Yl (2.4) 

7=1 

The complex eigenenergies and the expansion coefficients are determined by 

^ [H^-,\e) - E{e)N^^>]c^>{e) = o (7 = i, • • • ,7ms«), (2.5) 

7'=1 

where the overlapping and Hamiltonian matrix elements are respectively written as 

iV^y = (2.6) 

and 

H^^m = i^-ymOWy)- (2.7) 
The complex resonance energy is given, independently of 6 in principle, by 

Eres = Er — —, (2.8) 

where is the resonance energy with respect to the 3a breakup threshold and F 
is the total decay width. 

The symmetric three-body basis functions iZ'-y in (2.4) is written as 

= ^7(ri, Ri) + #^(r2, R2) + #7(r3, R3). (2.9) 

We express each <^7(rj, Rj) as a product of a function of rj and that of R^: 

^^{vi,ni) = ^nl{ri)i^NL{Ri)\Yi{vi)YL{%)] , (2.10) 

where 7 specifies a set of quantum numbers 

-i = {nl,NL,JM}. (2.11) 
J is the total angular momentum and M is its z-component. 
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2.2. Real-range Gaussian basis functions 

According to the Gaussian expansion method (GEM) 'US we take the 
radial shape of 4>ni{r) and ipNL{r) in (2.10) as follows: 

0nKO = ^'e-(^/'"^', (2.12) 
V^,VL(i?) = i?''e-(^/^^)', (2.13) 

where normalization constants are omitted for simplicity. The GEM recommends to 
set the Gaussian ranges in a geometric progression: 

r„ = ria"~^ (n = 1, • • • , n^ax) , (2.14) 
Rn = RiA''-^ (iV = l,--- ,iV^ax) . (2.15) 

This greatly reduces the nonlinear parameters to be optimized. We designate a set 
of the geometric sequence by {rimax, '^i, '^nmaxl instead of {n-max, fi, a} and similarly 
for {A'max; Ri, -^Afmax} ) which is more convenient for consideration of the spatial 
distribution of the basis set. 

The basis set {(pni'i n = 1, - ■ ■ , n-max} has the following properties: i) They range 
from very compact to very diffuse, more densely in the inner region than in the 
outer one. While the basis functions with small ranges are responsible for describing 
the short-range structure of the system, the basis with longest-range parameters is 
for the asymptotic behavior, ii) After multiplication by normalization constants for 
{4'ni I 4'ni) = 1) they have the relation 

/ 2a'' \ '^^^^ 

{(Pnl \ <Pn+kl) = L _^^2fc ) ' (2-^^) 

which shows that the overlap with the kth neighbor is independent of n, decreasing 
gradually with increasing k. We thus expect that the coupling among all the basis 
functions takes place smoothly and coherently so as to describe properly both the 
short-range structure and long-range decaying behavior simultaneously. 

We note that a single Gaussian e"^^/^") decays quickly as r increases, but appro- 
priate superposition of many Gaussians can decay even exponentially with increasing 
r up to a sufficiently large r. Good examples are shown in Figs. 3 and 4 in Ref. [9]). 

2.3. Complex-range Gaussian basis functions 

For the precise GSM calculations of three-body systems, however, we improve the 
Gaussian shape to have more sophisticated (but still tractable) radial dependence. 
This is because that the wave function in GSM becomes more oscillatory as the 
scaling angle 9 increases. But, the superposition of the real-range Gaussians is 
difficult to accurately describe oscillatory functions having several nodes. 

In the GEM of Ref. E]) , we proposed to improve the Gaussian shape by intro- 
ducing the complex range instead of the real one: 

^(+-)(^) ^^«g-(l+i-)(r/r„)2^ (2.17) 

0(^7^) (r) = r' e-(i-'^^)('-A")', (2.18) 
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and 

ij^+^\R) = i?^ e-^i+i-^^iR/RN)\ (2.19) 

V^^L ^(i?) = i?^ e-(i-^'^)(^/^^)', (2.20) 

where the ranges r„ and i?7v are given by (2.14) and (2.15), respectively. Using the 
above complex conjugate pairs, 4>^^^\r) and ip^^^\R), we can construct equivalent 
sets of real basis functions: 

^(cos)(^^ = e-("/''")' cosa;(r/r„)2, (2.21) 

</,Ji^)(r) = r' e-('^/'^")' sina;(r/r„)2, (2.22) 

and 

i^^^L\R) = i?^ e-^^/^'^)' cosojiR/RNf, (2.23) 



^^'°)(i?) = i?^ e-^^/^'^)' sinu;(i?/i?^)2. (2.24) 



/.(sin) 

■'nl 

In the present CSM calculation of ^^C, the former set, (2.17)-(2.20), is emplovecJ*^ 
and the three-body basis function (^^(rj,Ri) of (2.10) is replaced by 

^^(ri,R,) = </.if^(r)vSS:^(i?,) [Yi{?i}YL{Ri)]^^, (2.25) 
with 7 specifying a set 

7 = {±uj,nl; ±00, NL- JM}, (2.26) 

where one can take different w's between <j)(r) and ?/'(i?) although it was not necessary 
in the present 3a CSM calculation. 

The new basis functions (2.17)-(2.24) apparently extend the function space from 
the old ones, (2.12) and (2.13), since they have the oscillating components; their 
apphcations are seen, for example, in Refs.|9|,|19]),[20|),|2l])j22l)j23l)il0l). Note that the 
computation programming is almost the same as that for (2.12) and (2.13) although 
some of real variables are changed to complex ones. 

In order to compare visually the real-range and complex-range Gaussians, we 
plot, in Fig. 2, (Pni{r) of (2.12), (/'^^'^(r) of (2.21) and (/'^"^(r) of (2.22) with / = 0, 
r„ = 5 fm and a; = 1.0 and 7r/2. A good test of the use of complex-range Gaussians 
is to calculate the wave functions of highly excited states in a three-dimensional 
harmonic oscillator (HO) potential. We calculate the I = neutron wave function 
in the potential with huj = 15.0 MeV. The wave function, iP'o, is expanded in terms 
of totally 32 basis functions of the types (2.21) and (2.22): 



tmax 



Mr) = E + c^r^&\r)] . (2.27) 



ra=l 



The expansion coefficients and the eigenenergies are obtained by diagonalizing 
the Hamiltonian in the space. Nonlinear parameters of the complex-range Gaussian 



*-* We made the same calculation employing the latter set, (2.21)-(2.24), to crosscheck the 
computation programs and obtained, as a matter of course, the same result. 
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1.0 



- 0.0 



-0.5 



r e*<"'"' cos (fl(r/r„)^ 

r e""""''sin a(rlr„f 




r n= 5 fm 



10 

r[fm] 



15 



20 




Fig. 2. An example of the I — complex- 
range Gaussian basis functions presented 
in the form of Eqs. (2.21) and (2.22) with 
r-„ = 5 fm and uj — 1.0 and 7r/2. The 
functions are normalized to unity. 



Fig. 3. Wave function of the I = 0, 

36-quanta state for the HO potential using 32 
complex-range Gaussians (black dotted line) . 
Deviation from the exact one (red solid line) 
is everywhere less than 0.001 in the unit of the 
vertical axis. 



Table II. Energies (number of quanta) of the highly excited I — states of the HO potential 
calculated using 32 complex-range Gaussians with {ninax=^, ''i = 1.4fm, rn^^^ = 7.2 fm, oj = 
1.0}. Wave function of the 36-quanta state is illustrated in Fig. 3. 



Exact 


Complex-range 




Gaussians 





10-^^ 


12 


12.000000004 


16 


16.00000006 


20 


20.0000006 


24 


24.000005 


28 


28.00002 


32 


32.0003 


36 


36.0005 


40 


40.04 


44 


44.2 



set are {nmax = ri = 1.4 fm, rn,^^^^ = 7.2 fm, lo = 1.0}. In Table II, the 
calculated energies are compared, in a high accuracy, with the exact ones up to 
highly oscillating states. In Fig. 3, wave function of the 36-quanta state is compared 
with the exact one. The two curves for those wave functions overlap to each other 
everywhere; the difference is less than 0.001 in the unit of the vertical axis. 

We thus expect that use of the new basis set in three-body CSM calculations well 
describes the hightly oscillating wave functions of both the resonant and nonresonant 
continuum states even when the scaling angles becomes rather large. 
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§3. Application to 3a resonances in C 

3.1. Interaction of the 3a system 

We take the same model and interaction as those in Refs. and I12p . The 
potential Vaa is constructed by folding the effective A^-A^ interaction by Schmid- 
WildermuttP^ and the Coulomb potential into the density of the a cluster having 
the {Os)^ configuration. In Refs. [TT]) and [T2]) . the Vaa is adjusted to reproduce the 
experimental phase shift of the a-a system by taking l.OSxVaa- 

The Pauli principle between a clusters is taken into account by the OCM.'^^ 
The OCM projection operator)^ ^auU, in the Hamiltonian (2.1) is written by 

ypauii= hm AV|/)(/1, (3.1) 

A— s-oo ^ — ^ 
/ 

which rules out the Pauli-forbidden a-a relative states (/ = OS, IS, OD) from the 
three-body wave function. In this work, we take A = 10^ MeV. 

Since use of the 2a potential Vaa together with the Pauli potential T4>auii makes 
the energies of the ground-rotational-band states [0^, 2]^, 4^) lower than the observed 
values, the repulsive 3a potential V^a in (2-1) is introduced in Refs. [TT|) and [T2]) 
phenomenologically in the form 

V3airi,r2,rs) = V£ eM-l^{rl + rl + rl)], = 0.15 fm'^, (3.2) 

where V^^^ = 31.7 MeV, V^^ = 63.0 MeV and V^^ = 150.0 MeV are employed 
dependently on the total angular momentum J = 0^,2+ and 4+, respectively. 

3.2. The 0+ resonances 

Figure 4, taken from Ref. TT]), shows the O"*" eigenvalue distribution of the com- 
plex scaled Hamiltonian calculated with the real-range Gaussian basis functions. The 
scaling angle 9 = 16° was the maximum angle available in the calculation. The new 
broad 0^ state was predicted at Eres = 1-66 — iO.74 MeV, but this complex energy 
is not isolated from the [a + a + a]+pBe(0"'") + a] continuum states in Fig. 4 at 
6 = 16°. The energy was derived by an extrapolation based on the ACCC+CSM 
(see Fig. 2 of Ref. [rf])). 

Figure 5 illustrates the result of the present CSM calculation for the J = 0+ 
states at 6 = 16° and 26°. The low-energy part is enlarged in Fig. 6. All the nonlinear 
parameters used in the calculation are listed in Table III. The parameters for the 
Gaussian ranges are given in round numbers but further optimization of them does 
not significantly improve the present result; the same is for the ranges of the J = 2"*" 
and 4+ states. The slightly larger basis set for / = L = 4 is necessary to precisely 
generate the [^Be(4^) + a] continuum. 

The [a + a + a] continuum for 9 = 16° is much less scattered than that in 
Fig. 4 and the pBe(0"'') + a] continuum is distinguishable in Fig. 6 from the former 
continuum by ~ 0.1 MeV which corresponds to the energy of the '^Be(0''') resonance 
measured from the a-a threshold. It is to be emphasized that, in Fig. 5, the localized 
resonance poles of the 0^, 0^ and 0^ states for 9 = 26° remain at the same places 
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of them for 9 = 16° and that the complex energies of the three-body continuum at 
9 = 26° still form a nearly 'straight line' with little scattered members. The energies 
and widths of those resonances are almost the same as those obtained in Ref . Ilip . 




Re (E) [MeV] 

Fig. 4. The 0^ eigenvalue distribution of the complex scaled Hamiltonian for the 3q system 
obtained by Kurokawa and using the real-range Gaussian basis functions. The scaling 

angle is 9 — 16°. The state was predicted at Eres = 1.66 — iO.74 MeV, not localized from 
the [of + a + a] + [*Be(0^) + a] continuum (see the text). This figure is taken from Ref. 11) . 




Re (E) [MeV] 

Fig. 5. The 0''' eigenvalue distribution of 
the complex scaled Hamiltonian for the 3a 
system with the use of the complex-range 
Gaussian basis set in Table III. The scaling 
angles are 9 — 16° (black) and 26° (blue). 
This figure is to be compared with Fig. 4. 



o [a+a+a] continuum 




Re (E) [MeV] 



Fig. 6. The low-energy part of Fig. 5 is en- 
larged to show that the [*Be(0^) -I- q] con- 
tinuum is distinguishable from the [a-|-a-|-a] 
one by ~ 0.1 MeV which corresponds to the 
energy of *Be(0^) measured from the a-a 
threshold. 
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Table III. All the nonlinear parameters of the J = 0"*" three-body complex-range Gaussian basis 
functions for ^^C used in the GSM calculation of Fig. 5 at 6 = 16° and 26°. Total number of 
the basis is 7max = 3200 with uj = 7r/2. 
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40.0 
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30.0 
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40.0 


1152 



3.3. The new Ot resonance 



In order to investigate the new O'^ state that was predicted by Kurokawa and 
Katop^''^ we performed the CSM. calculation for the scaling angles from 9 = 22° 
up to 36° . These large angles are required to reveal explicitly such a low- lying broad 
resonance separated from the [a -|- a -|- a] -|- pBe(0"'") -|- a] continuum states. The 
employed set of the three-body complex-range Gaussian basis functions is listed in 
Table IV. Such a larger number of the basis is necessary for this purpose. 



0.0 



-0.5 



-1.0 



■Bo 




C (o^) 



9=22° 
/ 



22° 



26° 



0=36 



oOO ° 30° 



0.5 1.0 
Re (E) [MeV] 



1.5 



Fig. 7. The 0"*" eigenvalue distribution of the 
complex scaled Hamiltonian in which the 
angle 9 is varied from 22° to 36°. The Oj 
resonance appears, as the closed blue circle 
(36°), at Eres = 0.79 - iO.84 MeV. Only 
for 9 — 36° both the continua of [q + a + a] 
(open blue boxes) and [*Be(0^) -I- a] (open 
blue circles) are given, but the former is 
omitted for 9 < 36° for clarity of the figure. 
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Fig. 



8. Trajectory of the O3 state obtained 
by changing the strength parameter 5 of the 
auxiliary three-body potential (3.3). The 
blue closed circle at 5 = corresponds to 
the 0^ resonance in Fig. 7. The green box 
denotes the (ri_state predicted by Kurokawa 
and on the basis of the extrap- 

olation (the dash-dotted green curve, taken 
from Ref. HJl) using the AGGG+GSM. 



In Fig. 7, we illustrate the O"*" distribution of complex eigenvalues for 9 = 22° up 
to 36°. Only for 9 = 36° both the [a -|- a -|- a] continuum (open blue boxes) and the 
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Table IV. All the nonlinear parameters of the J = 0"*" three-body complex-range Gaussian basis 
functions for ^^C used in the GSM calculation of Fig. 7 for 6 — 22° to 36° to reveal the 
resonance state. Total number of basis is 7max = 4448 with u) — -k 12. 
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0.6 


40.0 
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22 
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40.0 


1936 


4 


12 
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30.0 
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12 
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['^Be(O^) +a\ continuum (open blue circles) are given, but the former is omitted for 
9 < 36° to avoid complexity of the figure. We observe a converged resonance pole at 
Eres = 0.79 — i 0.84 MeV (the closed blue circle) and identify it as the third 0"*" state 
that was predicted in Ref. [TT]) . The position and width of the resonance, however, 
differ slightly from the result of Ref. HU, Eres = 1-66 - i 0.74 MeV. 

Reason of this difference is explained with Fig. 8 which is to be compared with 
Fig. 2 of Ref. ,11). Figure 8 illustrates the trajectory of the 0'^ state on the complex 
energy plane, that was obtained by changing the strength parameter 5 of the auxiliary 
three-body potential, Eq. (4) of Ref. 11), added to the Hamiltonian (2.1), 

Kux. =<5exp[-;u(r2+r2 + ri)], ^ = 0.15 fm-^. (3.3) 

The closed blue circle for 5 = in Fig. 8 is the same as that for the 0^ resonance in 
Fig. 7. On the other hand, in Ref. llip . the direct CSM calculation of the 0^ resonance 
was not possible when the auxiliary 3a potential is less attractive than 6 = —120 
MeV. The green box that indicates the Oj^ state of Ref. [IT]) was therefore esti- 
mated by the extrapolation (the dash-dotted green curve) using the ACCC-I-CSM. 
We thus understand that the difference in the resonance-pole position between the 
two calculations comes from the error of the extrapolation. 

We conclude that we have confirmed the prediction by Kurokawa and KatcP^ 
about the appearance of a new Oj^ broad resonance that is located slightly above 
the Hoyle state (0^). As long as the structure of the Oj]" state is concerned, it 
is interesting to see that, in Fig. 7, the converged pole of the state is generated 
from the [^Be(0+) -|- a] continuum during the scaling angle is rotated up to = 36°. 
Therefore, the Oj^ state is considered to be dominantly composed of the [^Be(O^) -|-a] 
configuration. Kurokawa and KatcP^J pointed out that the 0^ state has a similar 
property to the 0^ state and the former may be a higher nodal state of the latter 
having the pBe(O^) -|- 0^(a)] configuration mainly, which is consistent with our 
conjecture. 

3.4. The 2"*" resonances 

Figure 9 illustrates the calculated 2^ eigenvalue distribution of the complex 
scaled Hamiltonian for the 3a system. The scaling angles are 9 = 16° and 26°. This 
figure is much more precise than Fig. 5 {9 = 16°) of Ref. [T?[) for the 2+ eigenvalue 
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Re (E) [MeV] 



Fig. 9. The 2 eigenvalue distribution of the 
complex scaled Hamiltonian calculated with 
the complex-range Gaussian basis set in Ta- 
ble V. The scaling angles are 6 — 16° (black) 
and 26° (blue). This figure is to be com- 
pared with Fig. 5 of Ref. [12]) at 61 = 16°. 
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Fig. 10. Low-energy part of the 2 eigenvalue 
distribution, where is varied from 26° to 
36°. Any 2'^ resonance, like the state in 
Fig. 7, does not appear in the upper-right 
side of the [a+a-t-a] continuum (blue boxes) 
ai 9 = 36° . See also the caption of Fig. 7. 



Table V. All the nonlinear parameters of the J = 2"*" three-body complex-range Gaussian basis 
functions for ^^C used in the GSM calculation of Fig. 9 at 6 — 16° and 26°. Total number of 
the basis is 7max = 6400 with u — -k 12. 
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distribution. All the nonlinear parameters used for calculating Fig. 9 are listed in 
Table V. Total number of the basis is 7max = 

6400 with uj = -k/2. 

Calculated five resonances denoted as 22,- ■■ ,2^ appear at almost the same 
complex energies of those obtained in Ref. I12|) . We observe no other resonance 
at low energies. As shown in Fig. 10, even if the scaling angle is increased up to 
9 = 36°, any new 2+ resonance, like the Oj^ resonance in Fig. 7, does not appear in 
the upper-right side of the [a -|- a -|- a] continumm at 6 = 36° . 
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3.5. The 4+ resonances 

The calculated 4"^ eigenvalue distribution of the complex scaled Hamiltonian 
is illustrated in Fig. 11. The resonance parameters are summarized in Table VI 
together with the result by Kurokawa and Kato.'^ All the nonlinear parameters 
of the basis set are listed in Table VII. We note that the lowest 4"*" resonance at 
Eres = 4.96 — i 1.1 MeV in Fig. 11 is missing in Ref. [T2|) . where the lowest one is 
given at Eres = 6.82— i 0.12 MeV, but it corresponds to the second 4+ state in Fig. 11. 




Re (E) [MeV] 

Fig. 11. The 4^ eigenvalue distribution of the complex scaled Hamiltonian for the Sa system 
calculated with the complex-range Gaussian basis set in Table VII. The scaling angles are 
9 — 16° (black) and 26° (blue). The lowest resonance at Eres ~ 4.96 — il.l MeV does not 
corresponds to the observed 4^ state at Er = 6.81 MeV which is reproduced by the second 4^ 
resonance in this figure. See the text about this problem. 



Table VI. Calculated resonance parameters for the 7 = 4"*" states in ^'^C by the present work and 
Ref. I12() which employ the same interaction. All quantities are given in MeV. 
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0.24 


4+ 


18.92 


11.62 
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17.11 
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This situation causes a serious problem in the determination of the 3a potential 
Vsa of Eq. (3.2). The strongly repulsive factor of V^^ = 150.0 MeV was so chosen 
in Ref. [T2|) that the calculated lowest 4+ state can reproduce the observed value of 
Eresi^i) ~ 6.808 — i 0.258 MeV. However, by the introduction of the repulsive 3a 
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Table VII. All the nonlinear parameters of the 7 = 4"*" three-body complex-range Gaussian basis 
functions for ^^C used in the CSM calculation of Fig. 11 at S = 16° and 26°. Total number of 
the basis is 7max ~ 8640 with uj = n/2. 



J 


= 4+ 




-{l±ii^)(r/r„)^ 






iui)(R/R 




ui = n/2 


I 


^max 


''I 


'■"max 


L 


max 


Ri 




number 






[fm] 


[fm] 






[fm] 


[fm] 


of basis 





20 


0.3 


30.0 


4 


18 


1.0 


40.0 


1440 


2 


20 


0.6 


30.0 


2 


18 


0.5 


40.0 


1440 


2 


20 


0.6 


30.0 


4 


18 


1.0 


40.0 


1440 


4 


20 


0.6 


30.0 





18 


1.5 


40.0 


1440 


4 


20 


1.0 


30.0 


2 


18 


1.0 


40.0 


1440 


4 


20 


1.0 


30.0 


4 


18 


1.5 


40.0 


1440 



potential, the original lowest 4"*" state at —1.66 MeV (with V^^ = 0) is pushed up 
to Eres = 6.61 — iO.lO MeV which corresponds to the observed 4^ state, but the 
second 4"'" state (with V^^ = 0) remains almost unaffected at Eresi^'^) = 4.96 — i 1.1 
MeV that becomes the lowest 4^ state in Fig. 11. 

Therefore, we understand that the introduction of the 3a potential V^a does not 
work for the J = 4"'" states even if the strength is given dependently on J. Any 
appropriate determination of the interaction in the 3a system will be required in 
future 3a OCM-CSM calculations. 

§4. Summary 

The authors have proposed to use the complex-range Gaussian basis functions, 
{r' e-i'^^^^)ir/rn) Yim(j)',rn in a geometric progression}, in the CSM calculations 
of many-body resonances in place of the real-range Gaussians that are often em- 
ployed in the literature. The former-type Gaussians are very suitable for describing 
short-range correlations, long-range asymptotic decaying amplitudes and highly os- 
cillating behavior in few-body systems as well as they are tractable in calculating the 
Hamiltonian matrix elements with transformation between different sets of Jacobi 
coordinates.!^ Therefore, they are particularly useful in the CSM calculations when 
representing the resonant and nonresonant continuum states that become quite os- 
cillatory as the scaling angle 6 increases; this enables us, in the study of broad 
three-body resonances, to take much larger angles than those considered before and 
to have a possibility of observing new broad resonance poles. 

The present method has been applied to the 3a resonances in ^^C with J = 
0^ , 2^ and 4"*" . The results were compared with those obtained by Kurokawa and 
KatdnJ'1^ where the real -range Gaussians were employed to expand the 3a wave 
function. In Table VIII, we summarize the calculated energies and widths of the 
states with J = 0"*" and 2"*" together with the result by Refs. [TT]) and [T2]) and the 
experimental data. The result for the 4+ resonances, having a problem in the inter- 
action employed, has been summarized in Table VI. 
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Table VIII. Summary of the calculated result for the J = 0"*" and 2'*' in ^^C by the present work 
together with the result by Refs.JJJ and il2)) and the experimental data. The data for Ex ~ 9.04 
MeV and 10.56 MeV are taken from Ref. 113^ : see the footnote in §1. The other data are taken 
from Ref. I26|) . The model and interaction are the same between the two calculations. All 
quantities are given in MeV. 
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The distribution of eigenvalues of the complex scaled Hamiltonian at ^ = 16° and 
26° were found to become very precise (Figs. 5, 6, 9 and 11). The prediction of the 
new broad Oj^ resonance by Kurokawa and KatcP^ was confirmed by our calculation; 
namely, as the scaling angle 9 increases up to 36°, the Ojj" resonance pole becomes 
clearly separated from the [a + a + a] and pBe(0"'") + a] continuum states (Fig. 7). 
The slight deviation of the Oj^ resonance energy by 0.9 MeV in Ref. .11,) from our 
result is attributed to the error of the extrapolatior^nj of the resonance position by 
the ACCC+CSM (Fig. 8). As for the 0|, 0^, 2+, 2+, 2+, 2+ and 2+ resonances, we 
obtained almost the same energies and widths as those in Refs. ITT|) and ll2|) . We did 
not observe any 2"*" low- lying broad resonance, like the Og", as long as we increased 
the scaling angle 9 up to 36° (Fig. 10). 

We employed the same interaction for the 3a system as used in Refs. [TT]1 and 
[T2]) . but the calculation did not satisfactorily well reproduce the observed energy 
of the important Hoyle state (0^), E^. = 0.38 MeV, with a deviation of some 0.4 
MeV higher. Furthermore, the strongly repulsive 3a potential for the 7 = 4"*" states, 
introduced in Ref. [T2|) to reproduce the observed energy of the 4^ state, is found 
to be not appropriate because the lowest 4"'" state in the present work is the broad 
resonance at Eyes = 4.96 — i 1.1 MeV (Fig. 11) which was missing in Ref. I12p and 
the second 4+ state corresponds to the observed 4^ state. 

We explicitly listed, in small tables, all the nonlinear parameters of the basis 
functions used in the present calculation of ^^C (0"'",2"'", 4"*"); our method is so 
transparent. For the comprehensive understanding of the 3a cluster structure of 
^^C, however, use of more improved interactions is highly desirable in future studies. 
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